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O ■ Abstract 
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. It has long been known that matter charged under a broken U(l) gauge symmetry collapsing to 

o ; 

form a black hole will radiate away the associated external (massive) gauge field. We show that the 
timescale for the radiation of the monopole component of the field will be on the order of the inverse 
Compton wavelength of the gauge boson (assuming natural units). Since the Compton wavelength 
for a massive gauge boson is directly related to the scale of symmetry breaking, the timescale for a 
black hole to lose its gauge field "hair" is determined only by this scale. The timescale for Hawking 
radiation, however, is set by the mass of the black hole. These different dependencies mean that 
for any (sub-Planckian) scale of symmetry breaking we can define a mass below which black holes 
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radiate quickly enough to discharge themselves via the Hawking process before the gauge field is 
radiated away. This has important implications for the extrapolation of classical black hole physics 
to Planck-scale virtual black holes. In particular, we comment on the implications for protecting 



protons from gravitationally mediated decay. 
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I. INTRODUCTION 



In the early 1970's Beckenstein and Teitelboim proved that a static bare black hole can 
be endowed with no exterior classical massive or massless scalar fields, nor with exterior 
classica, tnassive vector ne.ds Q. At the S a,„e thne, Price exp.ored the m echani Sm h y which 
matter collapsing to form a black hole could divest itself of a classical massless scalar field 
P). He concluded that the mass of the nascent black hole will set the timescale for a massless 
scalar monopole field to radiate away. 

This leaves open the question of the timescale for the loss of a massive vector field. The 
question is nontrivial, since setting the vector mass to zero results in the case of electromag- 
netism (unbroken U(l) gauge symmetry). It is well known that an external electric monopole 
field can persist even in the limit of a static black hole, where we recover the Reissner- 
Nordstrom solution. We are led, then, to consider whether a massive vector monopole will 
decay with a timescale set by the vector boson mass or whether (as in the case of a scalar 
field) the timescale is determined by the black hole mass. In the former case we recover 
electromagnetism in the continuous limit of small boson mass, while in the latter situation 
a massless photon represents a discontinuous jump from the physics of massive vector fields. 

Coleman, Preskill and Wilczek, in their work on discrete gauge symmetries and black 
holes, state (without proof) that the correct answer is the continuous one and that the 
lifetime of a classical massive vector monopole field is set by the mass of the field itself 3|. 
In this work, we show that their assertion is correct. We demonstrate this by adaptingthe 
detailed calculational techniques used by Price in his study of massless vector fields [3,0] to 
the case of a massive vector monopole field. 

This answer has important consequences for black hole phenomenology. Black holes are 
postulated to radiate via the Hawking process |5(. When charged particles fall into a black 
hole, this process can theoretically recover the "lost" charge if and only if some external field 
generated by those charges remains [fj, Q, 1^ ■ Since black holes cannot sustain an external 
massive vector field forever, it has been assumed that particles charged only under a broken 
U(l) gauge symmetry that fall into a black hole result in charge non-conservation - their 
charge will not be recovered through Hawking radiation. Now we see that this need not be 
the case for all black holes. The timescale for Hawking radiation is related only to the mass 
of the black hole. Therefore, because the external field's decay timescale and the black hole's 
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radiation timescale depend on parameters that are totally independent of one another, we 
can define a regime in which the black hole will discharge quickly enough to respect a broken 
symmetry. 

This result has particular relevance to the assumptions made about virtual (Planck-scale) 
black holes and the so-called "spacetime foam" (see, e.g. It has been assumed that since 
classical black holes do not respect a broken symmetry, quantum processes involving black 
holes will also violate such symmetries. In this work, however, we will show that the smallest 
classical black holes will in general respect broken symmetries. Thus, the extrapolation of 
conservation laws to quantum black holes must be re-examined. 



II. DESPUN FIELD EQUATION FOR A MASSIVE VECTOR BOSON 

We will consider the collapse of a star containingmatter that is charged under a broken 
U(l) symmetry by adapting the model defined in j^, 0] for the cases of stars acting as a 
source of massless scalar and vector fields to our purpose. We will be concerned only with 
the field external to the star (and later the black hole) so we will work in the Schwarzschild 
geometry with line element: 

ds 2 = (l- —) dt 2 -(l- —) 1 dr 2 - r 2 dd 2 - r 2 sin 2 Odcf 2 . (1) 



We will make the assumption that the charges in the collapsing star are small enough that 
the vector fields never have a significant impact on the external geometry. 

To analyze the behavior of the external massive vector field, it will be very convenient to 
use the "despun" versions of the field equations, as emphasized by Price [j| . This formalism 
is briefly reviewed in Appendix El I n the case of a purely monopole field, the field equations 
reduce to: 
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-d 2 r *% + d 2 %-d r * 
where r* is defined by: 

r* = r + 2Mln(^-l)+C (3) 

and C is an arbitrary constant. It is important to note that in these coordinates the horizon 
of the black hole will be found at r* = — oo. Note also that the quantity fi appearing in 
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Equation (J2J) is the inverse Compton wavelength of the vector boson field. We will assume 
throughout this paper that: 

liM < 1. (4) 

We will comment on this assumption in Section II VI For a monopole field, as explained 
in the Appendices, the despun quantity $o appearing in Equation (J2J) is simply the radial 
component of the electric field (or the equivalent for the massive vector boson field). As 
a check on the form of our field equation, then, it is important to notice that in the limit 
r ^> 2M Equation (J2J) admits the flat space Proca solution: 
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To make full use of Price's results, we must make one further simplification. We make 
the change of variables $o = which gives the convenient form: 
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III. CALCULATION OF DECAY TIMESCALE 



A. Reflection of the Field from the Curvature Potential 



We are now in a position to use the field equations of the previous section to determine 
how long a massive vector monopole field will persist outside of a collapsing star. In his 
analysis of scalar monopole fields, Price used both analytical and numerical methods to 
show that the timescale for decay is set by M, the mass of the collapsing star [2j. We will 
attempt to adapt both of his methods to the case of a massive vector field. 

We begin with the analytical approach. Price was able to show that the field equation 
for the despun piece of any radiatable massless field will reduce to the same form exhibited 
by Equation © - namely, the field equation will be a one-dimensional wave equation with 
an effective potential j^Q]. In the case of our massive vector monopole field, the effective 
potential is: 



eff 
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This effective potential plays an important role in the loss of the external field after 

flaur 



collapse of the star 



As seen in Figure ^ the barrier is highly localized at a 
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specific value of r* (corresponding to r ~ 3.5M). The precise value of r* depends on the 
choice of the arbitrary constant in Equation (J3j) . Waves propagating out from the surface of 
the star as it collapses beyond the location of the barrier will have to penetrate the potential 
barrier if the field far outside is to remain coupled to the sources within the star. 

As emphasized in Price's approach DDI 10j,lll(, we will be most interested in the behavior 



of long wavelength waves as they approach the barrier. This is because time dilation requires 
that the field on the stellar surface settle down to a constant value as: 

* = a + bex v( J -Jif) ® 

where a and b are constants. As the stellar surface approaches the horizon (t — > oo and 
r* — > — oo) the field becomes constant. This implies an effective redshift of the information 
propagating out from the star. If the external field is to remain in contact with the sources 
falling into the forming black hole, it can only do so via waves with wavelengths approaching 
infinity. We can find the transmission and reflection coefficients for such waves by examining 
the solution to the wave equation (0) in the three distinct regions (far outside the peak of 
the potential, the region where the potential is important, and well inside the peak of the 
potential) and then matching these solutions 10]. We will assume 

* = ip(r) exp(-iut). (9) 

We begin with the solution well outside the peak of the potential. In mathematical terms, 
we want the limit r* ^> 2M. Note that for very large values of r*, Equation Q implies 
r ~ r*. In this limit, the differential equation for ip(r) becomes: 
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^far + W^fcr 
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^far = (10) 



dr* 

which has a solution in terms of spherical Bessel functions of the first order To find 

the transmission coefficient for waves propagating out from the black hole, we must choose 
the linear combination which reduces to the form T exp(ikr*) in the limit of very large r* . 
This turns out to be: 

^ T =Th? (y^^r*y (ii) 

Next we look at the solution in the region where the potential is important. Recall that 
we are assuming /iM C 1 so that in this region the /i 2 contribution is negligible. Further, as 
we have stated, we are only interested in the reflection and transmission coefficients for long 



wavelength modes. Thus, upon assuming the time dependence of Equation we lose no 
information by taking the limit us <C M~ x . With all these simplifications, the field equation 
for ip near the peak of the potential becomes:: 



d 2 , / 2M 
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^peak = 0. (12) 
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This equation has exact solutions. One solution is easy to guess, as this equation is exactly 
the same field equation we would have found for a static, massless vector boson field. (The 
limits of small uj and small \i make this apparent.) This indicates that one valid solution is 
^peak ocl/r. The full solution near the peak can be written: 
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^peak = — + B 
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(13) 



where A and B are constants. 

Finally, we consider the region inside the potential. Here, we are in the limit of large 
negative values of r*. In this limit, the term (1 — 2M/r) causes the potential to vanish 
exponentially rapidly in r*. For this reason waves propagate freely near the horizon in r* 
coordinates. We can therefore adopt the solution (for outgoing waves): 

VWr = exp(ico>r*) + R exp(— iujr*) . (14) 

We now solve for the transmission and reflection coefficients by matching Equations (jllj) . 
(|13|) and (|14j) to form a complete solution. We begin by matching the far zone to the region 
near the potential peak. Recall that we are interested in the limit of long wavelength, so we 
can approximate the spherical Hankel function of Equation (fTTj) as: 

/ (^ 2 -^ 2 )(r*) 2 T T 

^far Z I . (15) 

where we have kept the leading order real and imaginary parts. Similarly, Equation (|13|) is 
of the form: 

^peak ~ 4 + ^ (16) 

r* 3 



Matching, we can see that: 



and: 



L= (17) 



B = T(lu 2 - /i 2 ). (18) 



Next, we match the near zone to the region near the peak. Here, for long wavelength, 
Equation (JT3jl becomes: 

i> neai ~l + R + icu(l-R)r*. (19) 



Equation (|13|) has the limit: 



Matching gives: 



and: 



VWk- ^ + 2MBr*. (20) 



A = 2M(1 + R) (21) 



B = ^-(l-R). (22) 
2M y ' y ' 

By equating our two expressions for A and our two expressions for B we can arrive at 
two equations in the unknowns T and R which are solved to obtain: 

2% 

T = } + ^ W^T) (23) 

and 

1 _ Ml! (w 2 _ 2)3/2 

l+4Mi (ci ;2_ 3) 8/2- ^ * 



We are left with a paradox. It seems that as u approaches u, the transmission coefficient 
drops to zero. Thus, we anticipate that, as discussed in j3,Q|, the external field will die 
away as the information from the source (the charge of the collapsing star) is reflected back 
and lost. It is important to note, however, that the transmission coefficient goes to zero in 
a way that is essentially independent of the precise value of fi (see Figure |2J). As long as fi 
is small compared to M, we see much the same behavior for T. This is a problem for the 
theory, since we find that even in the limit \x = taking u — > gives T = 0. This would 
mean that even in the case of unbroken U(l) symmetry the external gauge field can be lost. 

The loophole that saves the case of the photon (or other massless gauge boson) is con- 
tained in the setup of the problem rather than the solution. In the simultaneous limit \x = 
and io = 0, we do not actually have three regions worth of solutions to match to one another. 
In fact, when we take these limits, the exact solutions we derived for Equation (|13|) will be 
valid at any r. This is a very important point, as we shall see. 

The fact that the exception for massless bosons is not obvious leads us to be suspicious 
of the case of nonzero fi. The most dangerous case is uj — * fi, since here it seems that the 



potential in the far zone will strongly resemble that for the massless photon. There is an 
important difference, however. A nonzero mass requires that we retain the term 2M/z 2 /r in 
the potential. Thus, for uo ph fi, the equation for if) in the far zone becomes: 
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dr* 
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,2M 
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(25) 



This is in the form of a Coulomb wave equation, and can be solved with the Coulomb wave 
functions jl^ : 



ipiax = T (g^—tj, \/uj 2 - /i 2 r*) + iF^—i], \JuP- - /j? 



where 



/x 2 M 



(26) 



(27) 



\/uj 2 - fi 2 

When expanded near the interface with the region near the potential peak we find that this 
solution becomes (see Appendix [0 for a detailed derivation): 



iTiiM 1 ' 2 7r 1 / 2 e^/ 2 
3 ^/sinh(7r?y) 



far 



T ^smli^r?) 2 
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(28) 



This implies that we have replaced the A and B parameters of Equations (|17|) and (fTSj) with 
the expressions: 

T VsinhH) , /4 



and: 



£ = ^ /U M 1 / 2 4=^=(1 + r/Y 7 V - /i 2 ) 3 / 4 



(29) 



(30) 



a/ sinh(7rr]) 

Combining these expressions with Equations (|2"T|) and (j22j), which still hold, we can solve 
for T: 

"2/iM 3 / 2 Tr^W 2 



T 



and: 



^ ^/sinh(7rr/) 

1 a/ sinh(7rr/) 
2/iM 3 / 2 ^VV^ 

1 - i? 4/i 2 M 3 vre^ 



(1 + ^1/2^2 _ ^2)3/4 + 



n -1 



(l+7 7 2 )- 1 / 2 ( W 2 - /U 2 )- 1 / 4 



1 + i? c<j sinh(7rr/) 



(1 + r? 2 )^ 2 -^ 2 ). 



(31) 



(32) 
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These expressions for T and R show the distinction between a massless photon and a 
massive vector boson in their explicit \x dependence. If fi 7^ 0, then the limit of vanishing 
uj 2 — ji 2 gives: 

lim ( ^— - 1 = 8M 2 vru 5 M 5 . (33) 

This equation tells us that even for the case of vanishing uj 2 — [i 2 we can consistently define 
three regions of solution for ip as long as fi 7^ 0. This fact allows us to construct trans- 
mission and reflection coefficients for the barrier. Therefore, we can be confident that the 
transmission of massive vector waves really does approach zero as uj becomes small. There 
are no technicalities to save an external massive vector field. 

Now that we have decided that massive vector waves can be said to reflect from the 
curvature barrier, we must decide what that statement means. In particular, we need to 
understand the timescales for reflection and transmission. It is important to note that our 
construction of the transmission and reflection coefficients does not allow us to pinpoint the 
location where reflection takes place. Rather, we must consider the reflection to be a result 
of the accumulated effect of the potential barrier over its entire width. We can say with 
certainty, however, that it is the edge of the barrier that is farthest from the black hole which 
plays the decisive role. If we can consistently define such an edge as the place where the 
terms generated by the mass of the Proca field begin to dominate the curvature potential, 
we expect reflection and loss of the field. If we cannot define such an edge, then a static 
1/r 2 law holds and we expect to find a Reissner-Nordstrom black hole. 

We can estimate the approximate location of this edge by comparing the magnitude of the 
terms in the potential of the far zone. It is apparent from Equation (jUJ) that the changeover 
from a curvature dominated potential to a mass-dominated potential will happen when: 




which gives a characteristic location: 

rw// -1 . (35) 

To approximate the time taken for waves to begin reflecting, we can construct two separate 
arguments leading to the same basic result. The clearest argument is that unless u is smaller 
than fi, the /z term is negligible in the field equation. Thus, our first inclination that the 
field is not massless must come from waves with a characteristic timescale jU -1 . Therefore, 
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we expect the external field to persist (as in the massless case) until a time of order has 
elapsed. Another way to arrive at this conclusion is to consider what we have already said 
about the mechanics of reflection. Since it is the far edge of the potential at r ~ /x" 1 which 
will make the final difference between reflection and transmission, an outgoing wave cannot 
in effect be reflected until it has time to feel the effects of the entire potential barrier, all the 
way out to r ~ Assuming then that the wave travels near the speed of light (this will 
at least be good for a lower bound) the time that must elapse between the onset of collapse 
and the first reflection of outgoing waves originating near the collapsing surface is of order 

A*" 1 - 

B. Numerical Solution 

The previous section has motivated the idea that an external massive gauge field arising 
from sources that collapse into a black hole will decay with a characteristic timescale set by 
the mass of the gauge field. Now, we wish to provide a more concrete demonstration. To 
do so, we can numerically integrate the field equation for $ - 

The first step in this process is to establish the initial conditions for the integration. We 
will follow Price [2] in considering the case of a static star with initial radius Ri ~ 4M which 
suddenly begins to collapse into a black hole of mass M. This situation implies that at the 
moment when collapse begins, the field far away from the star should resemble the flat space 
static Proca solution. Thus, we expect that the "electric field" far from the star will have 
the form: 

(36) 

where Q is the total net Proca charge (in some normalized units) contained within the star. 
Since we know that $o is equal to the "electric field" for the Schwarzschild space Proca 
solution (see the discussion in Appendix El and around Equation ©), we can use the flat 
space form for E to set the field value and first derivative of $o & t some point far from the 
initial surface of the star. Then, by numerically integrating the field equation for $ outside 
the star (Equation ©), we can find the initial conditions for $ everywhere outside the star 
before collapse. 

To make this integration simpler, we will assume that our end result will resemble the 
flat space Proca solution everywhere (this assumption will be justified in our results). With 



E = Qe 
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this assumption, it becomes useful to make the substitution: 



$ 




(37) 



which results in the field equation for f(r): 




(38) 



In terms of f(r), the boundary conditions for our numerical integration are particularly 
simple. Far from the location of the star (r ^> M) we expect / = 1 and /' = 0. We choose 
to take /(5//i) = 1; /'(5//i) = as our far boundary condition and integrate in toward the 
initial radius of the star. The results of this integration are summarized in Figures 0HS1 

From these results, one can see that the initial conditions for $o just outside the star 
will deviate from the Proca solution at the order fiM. Now we must consider how the field 
evolves in the area vacated by the collapsing star once collapse to a black hole has begun 
(the region between Ri and the eventual horizon 2M). Here we must appeal to an analogy 
with electromagnetism. 

In the case of a zero mass field, $q M would be Q/r 2 everywhere outside the initial surface 
of the star. In this case, the Reissner-Nordstrom solution tells us that upon collapse of 
the star the $q M field will simply evolve to equal Q/r 2 everywhere external to the black 
hole. This is important, because numerical integrations for the Proca field show that if 
fi 2 R 2 ~ fi 2 M 2 << 1 we obtain initial conditions that closely resemble those for a massless 
EM field. In fact, our initial conditions on $o near the surface of the star can be summarized: 



where Q — Q (1 + O(fiRi)) renormalizes the EM profile (this is allowed because the overall 
normalization is irrelevant to the evolution of the EM field - we would obtain a Q/r 2 profile 
regardless of the value of Q). If we now take Ri of order 4M then we see from Figure Q 
that at radial distances r < Ri the effective potential governing the evolution of $o rapidly 
approaches that for electromagnetism (differing by less than order fi 2 M 2 ). Putting these 
initial conditions and the effective potential into the differential equation governing the 



d r $o(Q,Ri,ti) 
d 2 %(Q,R u ti) 



$o(Q,Ri,U) 



^ u (Q,Ri,ti) 

d r ®o M (Q, R4, U) {l + 0(fi 2 R 2 )) 

a r 2 $H4^) (i + o^ 2 r 2 )) 



(39) 
(40) 
(41) 
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evolution of the field (Equation (J2J)) shows us that the evolution of the <3>o held in the region 
vacated by the collapsing star will be governed by a differential equation that differs from 
that for a massless EM field only at the order fi 2 M 2 . Further, these differences will be 
suppressed by the factor (1 — 2M/r) which rapidly becomes small in this region. Similarly, 
the integration of the field equation will use boundary conditions that differ from the EM 
case only at the order fi 2 M 2 and the collapse will characteristically require a time M <C fi~ x . 
Taken together, all of this implies that we expect the profile of a massive vector monopole 
field to differ by 0(fi 2 M 2 ) from the EM case in the near neighborhood of the event horizon 
(r ~ M) immediately after the star has completed its collapse to a black hole. 

Our analogy with the EM case has shown us that the Proca field in the region vacated 
by the collapsing star will have a profile that is within order fi 2 M 2 of the EM field result 
Q/r 2 immediately following the collapse. This, in turn, differs only at the order fiM from 
the flat space Proca profile. We have already shown through explicit numerical integration 
that the field outside the original surface of the star will also be within order \xM of the flat 
space Proca profile. We can therefore assume that up to errors of order \xM <C 1: 



$„(f = 0) = Qexp {-fir) 



a 1 



(42) 



describes the Proca field profile everywhere outside the horizon of the black hole immediately 
after collapse. Since we are trying to show that the field will persist well after the black hole 
has collapsed, we will begin our simulations after the black hole has already formed. 

Using the static Proca field solution as our initial conditions, we numerically integrate 
the full field equation (J2J) for several values of fiM. The results of these integrations are 
summarized in Figures El and These figures clearly show that the timescale for loss of the 
field is characteristically of order and that the field near the black hole remains almost 
completely undisturbed for a time approximately equal to fi~ l as long as fiM < 1/40. For 
fiM > 1/40 we can make no definite statement about the field behavior, since in this regime 
the field evolution is sufficiently violent that a more sophisticated numerical integration 
technique is needed (also, one should begin to consider evolution of the field during the 
collapse of the star). 
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IV. CONSEQUENCES 



In this section we wish to discuss possible values of \iM. We therefore switch from the 
natural units of general relativity (G — h — c — 1) to the natural units of field theory 
{h = c = 1, G — 1/Mpt). The major difference is that M no longer has the units of radius. 
Instead, the Schwarzschild radius is 2M/Mp L . Thus, rather than the condition fiM <C 1, the 
correct expression in field theory units is: fiM/M^ <C 1. 

Once we have shown that a massive vector field around a star collapsing to a black hole 
will persist for a time of order /x -1 , we can draw an important conclusion. The lifetime of 

n 

a black hole goes as M 3 |5|. Since a Planck-mass black hole is expected to live one Planck 
time, we have the approximate relation: 

M 3 

tevap ~ Jjl ■ (43) 

This approximation is essentially borne out by the calculations of Page , up to a numerical 
factor which depends on the particle species that are emitted by the black hole. For the 
smallest black holes, one must include many known species not accounted for in Page's work 
and perhaps some wholly unknown species. These extra particles will tend to lessen the 
lifetime of small black holes, bringing the results of 3| closer to our own approximation. 

As long as t evap is less than /x" 1 , we anticipate that the external field will persist for the 
entire life of the black hole (unless, of course, the charge is recovered before the complete 
evaporation of the black hole). If the black hole can evaporate while the external field 
persists, we anticipate the recovery of the charge jf| Q, Isj]. We can therefore estimate the 
black hole mass below which a symmetry broken at the scale [i will obey charge conservation. 
The result is: 

M < -fc. (44) 

In a traditional universe with the Planck scale at about 10 19 GeV, this implies that a 
symmetry broken at a scale of 10 4 GeV (just beyond the reach of current colliders) would 
be respected by black holes up to 10 5 Planck masses (approximately 1 g). 
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V. EXTRAPOLATIONS 



It has historically been assumed that virtual or quantum black holes would not respect a 
broken symmetry because classical black holes do not. Here we have shown, however, that 
there is a strong case to be made for the claim that the smallest classical black holes will 
respect a broken symmetry. What constitutes a "small" black hole depends on the scale at 
which the symmetry is broken. Symmetries broken at small scales (high energies) require 
less massive (shorter-lived) black holes than symmetries broken at low energy. Since Planck- 
scale black holes are smaller and shorter-lived than any classical black hole, this result calls 
into question the conclusion that broken symmetries are not respected by quantum gravity. 

This result could be of particular importance if there is a low scale for quantum gravity. 

I I 

Here, virtual black holes have been predicted to result in fast proton decay [14]. Now, 
however, we see that protons might be protected by a broken gauge symmetry. In fact, it is 
conceivable that the scale for such symmetry breaking can be pushed nearly to the quantum 
gravity scale, leaving the possibility that it decouples from physics at currently explored 
energies even in a low Planck scale universe. 



APPENDIX A: DERIVATION OF THE FIELD EQUATIONS 

We begin our derivation of the field equations for a massive vector field in Schwarzschild 
spacetime by following Price in making use of the spinor formalism to write the field equation 
for a massive gauge boson in terms of scalar quantities. To obtain the correct expression 
'or the field tensor F^ v it is necessary to supplement the expressions given in Price's work 
"\ with the formalism found in the article by Pirani Q|. In this way we arrive at the 
expression: 

= a^a u B yF AkBt (Al) 

where Latin capitals denote spinor indices (see Q] for more detail). 

We can now choose the connection cr^ • to be determined by the null tetrad chosen to 

n 

parameterize our spacetime [lfij . We follow Price in using the null tetrad: 

F = ( (1-^) \l,0,0) (A2) 




n" = £ (1,-11- — 1,0,0) (A3) 



r 
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We then choose: 

<o = 1"; <i = n"l o& = ro"; er? = m» (A5) 
where mP 1 denotes the complex conjugate of m M . The conversion to scalar fields is now made 



via the definition |4|: 



Faxby = g OuB $ *y + c^ab) • (A6) 



Using the rules for contracting spinor indices (see [15J), we can express the field tensor 
in the form: 

F>» = ^[-P^($ 10 + $* ) + ZW($„i + + m"rr($ 10 - $Si) + cc 

+ (IV - rm^fu + c.c. + (m>W - m^n^oo + c.c] (A7) 

where c.c. denotes the complex conjugate of the term immediately preceding. Note that 
this expression for F^ v is manifestly antisymmetric and real as long as $10 = $01, which is 
guaranteed by the total symmetry of spin indices QQ. We desta an antisymmetric, real 
F^ v so that we can make the standard definition: 

= d^A v - d v A„. (A8) 

We henceforth simplify our notation by introducing: $ = $ 10 , $+1 = $oo> an d $-1 = $11 • 
We now have the tools necessary to derive the field equations for the Proca field. To 
obtain a convenient form, there is one more trick to use. We combine the homogeneous and 
inhomogeneous Proca equations (Maxwell equations for a massive vector field) into the form 

F^ v + -j=(e a ^F a ^), u = ^A". (A9) 
^ v 9 



This yields the equations (assuming $ has the form of a spherical harmonic): 

D (V 2 <l ) = r$ +1 + //V [A - S^^A 1 ] (A10) 

2A (V 2 l>o) = 2r&_i - rV [S{r)A° + A 1 ] (All) 

2D (r$-j) = -l(l + l)$ + fi 2 r 2 [i sin 6 A 3 - A 2 ] (A12) 

2A (rS(r)$ +1 ) = -1(1 + l)S(r)$ + (i 2 r 2 S{r) Q[A 2 + isinOA 3 } . (A13) 
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Here we have defined a number of quantities to simplify the form of the final equations. 
They are: 

D = (l-—) [d t + d r *] (A14) 



A = ^[d t - d r *\ (A15) 
6 = + (A16) 

Sir) ^ (l - . (A17) 

Further, we have employed the "despun" versions of the fields as defined in jj|: 

$o = $o (A18) 

$ +1 ee -L - -A-a + cot ^ $ +1 (A19) 

$_! EE -L ( 5 e + + C ot 0) $_l (A20) 

One benefit of using these definitions is that it is immediately obvious that the equations 
for a massive vector boson field reduce to Price's equations for a massless vector boson field 
(given as Equation (42) in 4]) when we take the limit /i — > 0. 



APPENDIX B: MONOPOLE FIELDS 

So far we have not assumed a monopole field. In fact, the results of Appendix [X] are 
general in the sense that any field can be decomposed into spherical harmonics. In this 
Appendix we will specialize to the case of spherical symmetry. The monopole produces the 
obvious simplification that I = 0. There are also restrictions on the form of F^ u imposed 
by spherical symmetry (assuming F^ u is given in terms of potentials as in equation (|A8|) ). 
Because the $ fields are defined in terms of F^ (see Equation ()A7jl ). spherical symmetry 
guarantees: 

$ = % (Bl) 
$ +1 = $_! = 0. (B2) 

Further, by comparing Equations (|A7|) and (|A8|) . we can see that in the spherically symmetric 
case: 

F 10 = $ . (B3) 
16 



By analogy with the case of electromagnetism, we see that $o parameterizes the effective 
"electric field" of the broken symmetry. 

Using these relationships, we can combine Equations (jAlOj) and (|A11J) into a single, 
second order differential equation. The resulting field equation is exactly Equation (J2J). 



APPENDIX C: COULOMB WAVE FUNCTIONS 



Equation (123)) is solved by the Coulomb wave functions Fi(— rj, p) and Gi(—r],p), where 
in our case 



V 



p 2 M 



and: 



p = v 1 OS 2 — p? r*. 



(CI) 



(C2) 



To consistently define a transmitted wave, we wish to find the linear combination of Fi 
and G_\ which asymptotes to e lkr * as r* (and hence p) approaches infinity. For p — > oo, we 
find (J: 

lim Fx(—r), p) = sin(p + (f>) (C3) 
where is a phase angle. Similarly, in the same limit: 



lim Gi(—r], p) = cos(p + (/>). 

p— >oo 

Thus, our transmitted wave will have the form: 

^ 6a = T(G 1 {-r } ,p)+iF 1 (-r } ,p)) 



(C4) 



(C5) 



just as asserted in Equation 

Now, we wish to look at the expansion of in the limit of small p, which will be 
matched to the solution for ip near the peak of the curvature potential (r* ~ few x M). For 
p <^ 1, we have 12]: 

F 1 (- v ,p)^C 1 (-r ] )p 2 (C6) 

where: 

exp(7rr / /2)|r(2-zr ? )| 



Ci(-V) = 

We can express |T(2 — irj)\ in the form 



12 



|r(2-^)| 2 = (i + ?/ 2 )- 



sinh(7r?7) 



(C7) 



(C8) 
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which gives: 



Similarly, for p <C 1, we can expand G\ to find [1 




V 



/2(i + 77 2)Vy. 



(C9) 



1 



(CIO) 



3Ci(-»?)p 



where Ci is the same function as above, giving: 




/2 (1+r/ 2 r i/: 



(Gil) 



These limiting forms give us exactly Equation (|28|). 
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FIG. 1: Effective potential of Equation Q as a function of r*, taking C = in Equation ©. The 
potential is plotted for three values of \x. The solid line represents \i = (unbroken symmetry), the 
dotted line [/, = 0.1, and the dashed line [i = 0.01. The line for n = 0.01 is almost indistinguishable 
from the massless case. Note that the plot for [i = 0.1 shows that the effects of \i on the potential 
are suppressed below r* = 3M (r ~ 3.5M). 
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FIG. 2: Magnitude of the transmission coefficient (\T\) for outgoing massive vector waves to 
propagate through the effective potential barrier as a function of the wave frequency u). Plots 
for the values fiM = 0.01, 0.001, and 0.0001 are overlaid. Note that the curves are essentially 
indistinguishable. 
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FIG. 4: Results of integration for f'(r) of Equation I381 This figure shows the results for fj, = 0. 
and // = 0.001 overlaid (indistinguishable). 
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FIG. 5: Results of integration for f"(r) of Equation I38I This figure shows the results for \i = 0.01 
and fM = 0.001 overlaid (indistinguishable). 
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FIG. 6: Evolution of the $ field for various values of fiM: 0.001 (+); 0.0025 (*); 0.005 (•); 0.0075 
(0); 0.01 (A); 0.025 (□). 
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FIG. 7: Log plot of the evolution of $ - for various values of juM: 0.001 (+); 0.0025 (*); 
0.005 (•); 0.0075 (0); 0.01 (A); 0.025 (□). 
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